The toroidal momentum pinch velocity 
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In this letter a pinch velocity of toroidal momentum is shown to exist for the first time. Using 
the gyro-kinetic equations in the frame moving with the equilibrium toroidal velocity, it is shown 
that the physics effect can be elegantly formulated through the "Coriolis" drift. A fluid model is 
used to highlight the main coupling mechanisms between the density and temperature perturbations 
on the one hand and the perturbed parallel flow on the other. Gyro-kinetic calculations are used 
to accurately asses the magnitude of the pinch. The pinch velocity leads to a radial gradient of 
the toroidal velocity profile even in the absence of a torque on the plasma. It is shown to be 
sizeable in the plasmas of the International Thermonuclear Experimental Reactor (ITER) leading 
to a moderately peaked rotation profile. Finally, the pinch also affects the interpretation of current 
experiments. 
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In a tokamak the total toroidal angular momentum 
is a conserved quantity in the absence of an external 
source. Transport phenomena determine the rotation 
profile which is of interest because a radial gradient in 
the toroidal rotation is connected with an ExB shearing 
that can stabilise micro-instabilities [3, 0, Q and, hence, 
improve confinement. Furthermore, a toroidal rotation of 
sufficient magnitude can stabilise the resistive wall mode 
0) IE El ■ In present day experiments the rotation is often 
determined by the toroidal torque on the plasma that re- 
sults from the neutral beam heating. Such a torque will 
be largely absent in a reactor and it is generally assumed 
that the rotation, and hence its positive influence, will be 
small. The novel pinch velocity described in this letter, 
however, may generate a sizeable toroidal velocity gra- 
dient in the confinement region even in the absence of a 
torque. 

We will focus on the Ion Temperature Gradient (ITG) 
mode, which is expected to be the dominant instability 
governing the ion heat channel in a reactor plasma. The 
equations are formulated using the gyro-kinetic frame- 
work 0, H, H which has been proven success- 
ful in explaining many observed transport phenomena 
Q|lll^|ll^|ia|lllllllll^|2E|2|. Because 
of the rotation, the back grou nd electric field cannot be 
ordered small [H, EH, [25L |26|] , and the starting point for 
the derivation is a set of equations for the time evolution 
of the guiding centre X and the parallel (to the magnetic 
field) velocity component (vy) in the co-moving system 
(with background velocity uq) obtained from Ref. [26[ 



dX 
~d7 



d«j| 
di 



b 

pR* 

B* 



mB 



x (eV^ + ^VB + mUo • Vuq), (1) 



(e V> + pSJB + muo • Vuo ) . (2) 



Here b = B/£? is the unit vector in the direction of the 
magnetic field (B), 4> is the perturbed gyro-averaged po- 
tential (i.e. the part not connected with the background 



rotation), p, the magnetic moment, m (e) the particle 
mass (charge), and Ug = Uq + vyb. For the background 
velocity (uo) we assume a constant rigid body toroidal 
rotation with angular frequency f2 ( this is an equilibrium 
solution see, for instance, Refs. [la US El] ) 



u = ft x X = R 2 nv<p, 
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where ip is the toroidal angle. We briefly outline the 
derivation of the final equations here. More details can 
be found in [29| . The background velocity Uo will be as- 
sumed smaller than the thermal velocity, and only the 
terms linear in Uo will be retained. This eliminates the 
centrifugal forces but retains the Coriolis force. Further- 
more the low beta approximation is used for the equi- 
librium magnetic field (i.e. b • Vb w V±B/B where _L 
indicates the component perpendicular to the magnetic 
field). With these assumptions 



B 



2v«n x b. 



(4) 



Using the definition of B* (see Ref. |26[) and expanding 
up to first order in the normalised Larmor radius p* = 
p/R, where R is the major radius, one obtains 
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and Sjf = b • B* = B(l + 2fi||/u; c ) {u) c = eB/m is the 

gyro- frequency) . Expanding now the equations of motion 
retaining only terms up to first order in p* yields 
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The terms in this equation are from left to right, the par- 
allel motion (viib), the ExB velocity ve, the combination 
of curvature and grad-B drift v^, and an additional term 
proportional to An interpretation of this term can 
be found if one uses the standard expression for a drift 
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velocity (vd) due to a force (F) perpendicular to the 
magnetic field vp = F x B/e£? 2 . Substituting the Cori- 
olis force F c = 2mv x ft, and taking for the velocity (v) 
the lowest order (parallel) velocity one obtains 
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The last term in Eq. ([6]) is therefore the Coriolis drift. 
Expanding the terms in the equation for the parallel ve- 
locity to first order in p* one can derive 
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where dX/di is given by Eq. The derived equations 
are similar to the non-rotating system, with the differ- 
ence being the additional Coriolis drift. It follows that 
this Coriolis drift appears in a completely symmetric way 
compared with the curvature and grad-B drift. 

In this letter the approximation that assumes circular 
surfaces and small inverse aspect ratio (e) is used. In this 
case the Coriolis drift adds to the curvature and grad-B 
drift 
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where e z is in the direction of the symmetry axis of the 
tokamak. The linear gyro- kinetic equation is solved us- 
ing the ballooning transform [30]. The equations, except 
from the Coriolis drift are standard and can be found in, 
for instance, Ref. [3l|. In the following vl = —KSJRVt/vth 
and u = RCl/vth- Unless explicitly stated otherwise all 
quantities will be made dimcnsionlcss using the major ra- 
dius R, the thermal velocity v t h = y/2T/rrii, and the ion 
mass rrii. Densities will be normalised with the electron 
density. The toroidal momentum flux is approximated by 
the flux of parallel momentum (T^) which is sometimes 
normalised with the total ion heat flow (Qi) 
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where / is the (fluctuating) distribution function and the 
brackets denote the flux surface average. 

Before turning to the gyro-kinetic calculations, first 
the implications of the Coriolis drift will be investigated 
using a simple fluid model (more extended models have 
been published in Refs. [12, HH). A (low field side) slab 
like geometry will be assumed with all plasma parameters 
being a function of the x-coordinate only. The magnetic 
field is B = Be y , VB = —B/Re x , The model can be 
build by taking moments of the gyro-kinetic equation in 
(X, uii , v±) coordinates 
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where Fm is the Maxwell distribution. Note that trans- 
lation symmetry in the z-direction is assumed, eliminat- 
ing the parallel dynamics. Building moments of these 
equations neglecting the heat fluxes (this a clear simpli- 
fication, see for instance [H, [H, [H, [37j ) , and taking the 
space and time dependence of the perturbed quantities 
as exp[ik$z — iuit], one arrives at the following equations 
for the perturbed density (n) normalised to the back- 
ground density (no), the perturbed parallel velocity (w) 
normalised with the thermal velocity, and the perturbed 
ion temperature (T) normalised with the background ion 
temperature (To) 
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Here R/L N = -i?Vn /n , R/L T = -RVT /T , the po- 
tential 4> is normalised to Tb/e, and the frequency is nor- 
malised with the drift frequency lud = —koTo/eBR. The 
Coriolis drift (all terms proportional to u) introduces the 
perturbed velocity in the equations for the perturbed 
density, and temperature. However, since u <C 1 the 
influence of the Coriolis drift on the "pure" ITG (with 
u = 0) is relatively small. The Coriolis drift generates 
a coupling between w and the density, temperature as 
well as potential fluctuations. Note that for u = the 
perturbed velocity is directly related to the gradient it', 
resulting in a purely diffusive flux. For finite rotation 
(u ^ 0) the ITG will generate a perturbed parallel veloc- 
ity w, which is then transported by the perturbed ExB 
velocity. If the perturbed temperature is kept the expres- 
sions for the momentum flux become rather lengthy and 
are, therefore, reported elsewhere [2^|. Retaining only 
the coupling with the perturbed density and potential, 
and assuming an adiabatic electron response (n = 4>/t 
with r = T e /Tq being the electron to ion temperature 
ratio) one can derive 
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Here, the dagger denotes the complex conjugate, lur is 
the real part of the frequency, and 7 the growth rate 
of the mode. Note that x<t> ' IS positive since ojr (7) are 
normalised to uj£> = —kTo/eBR. The second term be- 
tween the square brackets of Eq. (|15[) represents an in- 
ward pinch of the toroidal velocity (the word pinch is 
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used here because the flux is proportional to u, unlike 
off-diagonal contributions that are due to pressure and 
temperature gradients [38|,[39j]) If one assume no torque, 
i.e. = it can be seen that the pinch can lead to a 
sizeable gradient length R/L u = RVu/u = 4 (for r = 1). 
The peaking is in roughly the same range as the expected 
density peaking |40| ]. 




FIG. 1: {R/2LT)T^/Qi as a function of u for three values of 
kepi 0.5 (o), 0.2 (squares), and 0.8 (diamonds). The top right 
graph shows the growth rate as a function of u and the down 
left graph the contour lines of (R/2LT)T$/Qi as a function 
of u and u' , both for kepi = 0.5. In the latter graph the thick 
line denotes zero momentum flux, i.e. the stationary point 
for zero torque 

Fig. [1] shows the parallel momentum flux as a func- 
tion of the toroidal velocity u obtained from line ar g yro- 
kinetic calculations using the LINART code [41| (in 
which unlike Eq. (| 1 1 [) the parallel dynamics is kept) 
for three different values of the poloidal wave vector 
(kg pi — 0.2, 0.5, and 0.8). The parameters of each of 
the gyro-kinetic calculations in this letter are those of the 
Waltz standard case [i^ : q — 2, magnetic shear s = 1, 
e = 0.1, R/L N = 3, R/L T = 9, r = 1, u = u' = 0. 
In the presented scans one of these parameters is var- 
ied while keeping the others fixed. Since the flux from 
Fig. [1] is linear in the velocity, a constant pinch velocity 
exists in agreement with the fluid model. The influence 
of the toroidal velocity on the growth rate is small. The 
bottom left graph shows the contour lines of the flux as 
a function of u and vl . The fact that the contour lines 
are straight means that the momentum flux is a linear 
combination of the diffusive part (cx x<P u ') an d t ne pinch 
velocity (V^u) 
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The diag onal part has been calculated previously using 
flui d lia 0, [HI, HI, S3 as weu as gyro-kinetic theory 
[48l l49j . The pinch velocity is negative (inward) for posi- 
tive u such that it enhances the gradient. It changes sign 
with u such that for negative velocities it will make u' 



more negative, i.e. the pinch always enhances the abso- 
lute value of the velocity gradient in agreement with the 
results from the fluid theory. Fig. Q] also shows that the 
pinch decreases with kgpi. It is noted here that also \<t> 
in becomes smaller for smaller kgpi [39j. 




FIG. 2: RV$lx<t> as a function of various parameters: R/Ln 
(x), 3s (+), q (o), and 20e (diamonds), and Qkepi (squares) 

Fig.[2]shows the normalised pinch velocity RV^jx^ as a 
function of various parameters. The magnetic shear and 
the density gradient have a rather large impact. Note 
that both due to s, as well as due to g, R/L^ and e, the 
pinch velocity is expected to be small in the inner core, 
but sizeable in the confinement region. 

The novel pinch velocity described in this letter has 
several important consequences. It can explain a gradi- 
ent of the toroidal velocity in the confinement region of 
the plasma without momentum input. A spin up of the 
plasma column without torque has indeed been observed 
[HI [H [H [H HE ill. Although & consistent descrip- 
tion ordering the different observations is still lacking, the 
calculations of this letter show that the pinch velocity is 
expected to play an important role. This finite gradient 
without torque is especially important for a tokamak re- 
actor in which the torque will be relatively small. From 
the calculations shown above, and for typical parameters 
in the confinement region of a reactor plasma, one ob- 
tains a gradient length R/L u — u'/u in the range 2-4 
representing a moderate peaking of the toroidal velocity 
profile similar to that of the density. Unfortunately, the 
current calculation only yields the normalised toroidal 
velocity gradient. In order to determine the velocity gra- 
dient one would need to know the edge rotation velocity. 
This situation is similar to that of the ion temperature 

The existence of a pinch can resolve the discrepancy 
between the calculated x<t> an d the experimentally ob- 
tained effective diffusivity (xcfi — T^/u'). The latter is 
often found to decrease with increasing minor radius and 
to be smaller than the theoretical value of x<t> in the outer 
region of the plasma HI, m • The pinch indeed leeds 
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to a decrease of \eS 
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The calculations in this letter show that the second term 
in the brackets can be of the order -1, leading to Xctf < Xi- 



(18) 



[32] 



R.E. Waltz et., Phys. Plasmas 1, 2229 (1994) 

T.S. Hahm, K.H. Burrell, Phys. Plasmas 2 1648 (1995) 

H. Biglary et al., Phys. Fluids B 2 1 (1990) 

A. Bondeson et al., Phys. Rev. Lett. 72, 2709 (1994). 

E.J. Strait et al., Phys. Rev. Lett. 74, 2483 (1995). 

H. Reimerdes et al., Phys. Plasmas 13, 056107 (2006) 

E.A. Friedman et al., Phys. Fluids 25, 502 (1982) 

D. H.E. Dubin et al., Phys. Fluids 26 3524 (1983) 
W.W. Lee, J. Comput. Phys. 72 243 (1987) 
T.S. Hahm, Phys. Fluids 31, 2670 (1988) 

C. Bourdelle et al., Nucl. Fusion 42, 892 (2002) 

E. J. Synakowski et al., Plasma Phys. Contr. Fusion 44 
A165 (2002) 

J. Candy et al., Phys. Rev. Lett. 91, 045001 (2003) 
X. Garbet et al., Nucl. Fusion 43, 975 (2003) 

D. R. Ernst et al., Phys. Plasmas 11 2637 (2004) 
M. Romanelli et al., Phys. Plasmas 11 3845 (2004) 
J.E. Kinsey et al., Nucl. Fusion 45, 450 (2005) 
A.G. Peeters et al., Phys. Plasmas 12 022505 (2005) 

F. Jenko et al., Plasma Phys. Contr. Fusion 47, B195 
(2005) 

C. Angioni et al., Phys. Plasmas 12, 112310 (2005) 
A.G. Peeters et al., Nucl. Fusion 45, 1140 (2005) 
A.Bottino et al., Plasma Phys. Contr. Fusion 48, 
(2006) 

T.S. Hahm, Phys. Fluids B4, 2801 (1992) 
M. Artun, Phys. Plasmas 1, 2682 (1994) 
T.S. Hahm, Phys. Plasmas 3, 4658 (1996) 
A.J. Brizard, Phys. Plasmas 2, 459 (1995) 
S.P. Hirshman et al., Nucl. Fusion 21, 1079 (1981). 
A.G. Peeters, Phys. Plasmas 5, 763 (1998) 
A.G. Peeters. et al., The toroidal momentum pinch 
be submitted to Phys. Plasmas (2007) 
J.W. Connor et al., Phys. Rev. Lett. 40, 396 (1978). 
M. Kotschenreuther et al., Comput. Phys. Commun 
128 (1995). 

J. Weiland, et al., Nucl. Fusion 29, 1810 (1989) 



215 



to 



88, 



[33] R.E. Waltz et al., Phys. Plasmas 4 2482 (1997) 
[34] W. Dorland et al., Phys. Fluids B 5, 812 (1993) 
[35] G.W. Hammett et al., Plasma Physics Contr. Fusion 35, 
973 (1993) 

[36] M.A. Beer et al., Phys. Plasmas 3, 4046 (1996) 
[37] Bruce D. Scott, Phys. Plasmas 12, 102307 (2005) 
[38] B. Coppi, Nucl. Fusion 42, 1 (2002) 

[39] A.G. Peeters et al., Plasma Phys. Contr. Fusion 48, B413 
(2006). 

[40] C. Angioni et al., Phys. Rev. Lett. 90, 205003 (2003) 
[41] A.G. Peeters et al., Phys. Plasmas 11 3748 (2004) 
[42] R.E. Waltz et al., Phys. Plasmas 2, 2409 (1995) 
[43] N. Mattor et al., Phys. Fluids 31, 1181 (1988) 
[44] S.-I Itoh, Phys. Fluid B 4, 796 (1992) 
[45] R.R. Domingruez et al., Phys Fluid B 5, 3876 (1993) 
[46] P.H. Diamond et al, Proceedings of the 15th IAEA Con- 
ference on Plasma Physics and Controlled Nuclear Fusion 
Research. (Sevilla 1994) p. 323 (IAEA Vienna 1994) 
[47] X. Garbet et al., Phys. Plasmas 9, 3893 (2002) 
[48] A.G. Peeters et al., Phys. Plasmas 12, 072515 (2005). 
[49] J.E. Kinsey et al., Phys. Plasmas 12, 062302 (2005) 
[50] M. Kotschenreuther et al., Phys. Plasmas 2, 2381 (1995) 
[51] L.G. Eriksson et al., Plasma Phys. Contr. Fusion 39, 27 
(1997) 

[52] J.E. Rice et al., Nucl. Fusion 39 1175 (1999) 
[53] I.H. Hutchinson et al., Phys. Rev. Lett. 84, 3330 (2000) 
[54] J.E. Rice et al., Nucl. Fusion 44 370 (2004) 
[55] J.S. deGrassie et al., Phys. Plasmas 11, 4323 (2004) 
[56] A. Scarabosio et al., Plasma Phys. Contr. Fusion 48, 663 
(2006) 

[57] D. Nishijima et al., Plasma Phys. Contr. Fusion 47, 89 
(2005) 

[58] P.C. de Vries et al., Plasma Phys. Control. Fusion 48, 
1693 (2006) 

[59] C. Angioni et al., "Theoretical understanding of core 
transport phenomena in ASDEX Upgrade", to be sub- 
mitted to Nucl. Fusion. 



